Two-dimensional surface charge transport in topological insulators 
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We construct a theory of charge transport by the surface states of topological insulators in three 
dimensions. The focus is on the experimentally relevant case when the Fermi energy ef and trans- 
port scattering time r satisfy eft jYi ^ 1, but ef lies below the bottom of the conduction band. 
Our theory is based on the spin density matrix and takes the quantum Liouville equation as its 
starting point. The scattering term is determined to linear order in the impurity density Ui and 
explicitly accounts for the absence of backscattering, while screening is included in the random phase 
approximation. The main contribution to the conductivity is oc n~^ and has different carrier density 
dependencies for different forms of scattering, while an additional contribution is independent of n^. 
The dominant scattering angles can be inferred by studying the ratio of the transport time to the 
Bloch lifetime as a function of the Wigner-Seitz radius rs- The current generates a spin polarization 
that could provide a smoking-gun signature of surface state transport. We also discuss the effect on 
the surface states of adding metallic contacts. 
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I. INTRODUCTION 

Topological insulators are a new class of materials that 
have been intensely researched in recent years.'!' They are 
band insulators in the bulk while conducting along the 
surfaces, possessing surface states with a spin texture 
protected by time-reversal symmetry. Within a certain 
parameter range the surface states are well described by 
a Dirac cone, allowing for parallels with graphene and 
relativistic physics, and prohibiting backscattering. Cer- 
tain physical properties related to the surface states of 
these materials can be expressed in terms of topologi- 
cally invariant quantities, and are thus protected against 
smooth perturbations invariant under time-reversal. For 
this reason topological insulators are being explored with 
a view towards applications, as a potential platform for 
quantum computation,'^' and as rich physical systems in 
their own right. 

The concept of a topological insulator dates back to the 
work of Kane and Mele, who focused on two-dimensional 
systems."^ These authors demonstrated that, in addition 
to the well-known TKNN invariant,^ characterizing the 
time-reversal breaking quantum Hall state, an additional 
topological invariant v exists, which characterizes the 
surface states of time-reversal invariant insulators. In 
two dimensions this invariant can be zero or one, the for- 
mer representing ordinary insulators (equivalent to the 
vacuum) and the latter topological insulators. The two 
band structures corresponding to i^ = and v = \ cannot 
be deformed into one another. The invariant v, which is 
related to the Z2 topology, counts the number of pairs 
of Dirac cones on the surfaces. Kane and Mele^ pro- 
posed realizing a topological insulator in graphene by 
using the spin-orbit interaction to open a gap, yet the 
spin-orbit interaction in graphene is rather small. Sub- 
sequently Bernevig et al. predicted that a topologically 
insulating state, the quantum spin Hall insulator, could 
be realized in HgCdTe quantum wells. ^ This state was 



observed shortly after its prediction by Konig et alW (for 
a review of this effect see Ref. lZj) The quantum spin 
Hall state constitutes a topological insulator in two di- 
mensions. Soon afterwards topological insulators were 
predicted to exist in three dimensions ,^ ^^ ^^ in which 
four Z2 topological invariants exist, customarily labeled 
I'D I i^ij 1^2 and 1^3. The Z2 invariants can be calculated 
easily if the system has inversion symmetry.lSJ The invari- 
ant vq determines whether a topological insulator is weak 
(i^o = 0) or strong {vq = 1.) The absence of backscatter- 
ing gives rise to phenomena such as Klein tunneling and 
is believed to make strong topological insulators immune 
to Anderson localization.^ 

A significant fraction of theoretical research on these 
materials to date has focused on their topology. A num- 
ber of articles have been devoted to the classification of 
topological insulators.'i2"21 A description of topological 
insulators in terms of topological field theory has been 
given in Ref. 15. and a general theory of topological in- 
sulators and superconductors has been formulated in Ref. 
[TCI An effective continuous model for the surface states 
of three dimensional topological insulators has been pro- 
posed in Ref. |T71 Their effective action in the presence of 
an electromagnetic field has been studied with a focus 
on the electric-magnetic dualityjE^ and analogies have 
been made with the topological invariants of the Stan- 
dard Model of particle physics-li^ The surface states of 
topological insulators in strong magnetic fields have been 
investigated in Ref. l20l 

The interface between a topological insulator and 
a superconductor was predicted to support Majorana 
fermions,'^ that is, particles that are their own an- 
tiparticles. Subsequent articles have focused on prob- 
ing Majorana fermions,'22l including a scheme employing 
interferometry.'^-J Recently the proximity effect between 
a superconductor and a topological insulator has been 
re-examined .ISS In addition, other exotic states have also 
been predicted, including a helical metabl^S supercon- 



ductivity and ferromagnetism induced by the proximity 
effect }2£l unconventional superconductivityj^ a quantized 
anomalous Hall effect in topological insulators doped 
with magnetic impurities,^ and a new critical state 
attributed to the Coulomb interaction.^^ Further ex- 
pected effects include the possibility of Cooper pair in- 
jection at the interface between a topological insulator 
and a superconductor j221 the inverse spin-galvanic effect 
at the interface between a topological insulator and a 
ferromagnet p^ giant magnetoresistance™ and a giant 
Kerr effect P 

Variations have also appeared on the theme of topo- 
logical insulators. Given that superconductors are also 
gapped, topological superconductors have been pre- 
dicted to exist .'i'^^iSlEl] Developments in different direc- 
tions include the topological Anderson insulator ,E2'33 the 
topological Mott insulator,'^ the topological magn etic 
insulator,^" the fractional topological insulator,SiIl2l the 
topological Kondo insulator,"^ as well as a host of topo- 
logical phases arising from quadratic, rather than linear, 
band crossings.'^ Topological phases in transition-metal 
based strongly correlated systems have been the subject 
of Refs. 1451461 and a theory relating topological insula- 
tors to Mott physics has been proposed in Ref. HT] 

Following the discovery of the quantum spin-Hall 
state, several materials were predicted to be topologi- 
cal insulators i n thr ee dimensions. The first was the 
alloy Bii_^Sb:„I^Sll9l followed by BizSeg, BiaTeg, and 
Sb2Te3.'^In particular Bi2Se3 and Bi2Te3 have long been 
known from thermoelectric transport as displaying siz- 
able Peltier and Seebeck effects, and their high quality 
has ensured their place at the forefront of experimental 
attention.'!' Initial predictions of the existence of chiral 
surface states were confirmed by first principles studies of 
Bi2Se3, Bi2Te3, and Sb2Te3.E3 On an interesting related 
note, it has recently been proposed that topologically in- 
sulating states can be realized in cold atoms,!^ though 
these lie beyond the scope of the present work. 

After a timid start, experimental progress on the ma- 
terials above has skyrocketed. Hsieh et al^using angle- 
resolved photoemission spectroscopy (ARPES), were the 
first to observe the Dirac cone-like surface states of 
Bii-jjSb^:, with x w 0.1. Following that, the same 
group used spin-resolved ARPES to measure the spin 
polarization of the surface states, demonstrating the cor- 
relation between spin and momentum.'^ Shortly after- 
wards exper iments identified the Dirac cones of Bi2Se3,i^ 
Bi2Te3,'S^157] a,nd Sb2Te3.E3 Recently Bi2Se3 nanowires 
and nanoribbons have also been growrpS and experi- 
ments have begun to investigate the transition between 
the three- and two-dimensional phases in this material.ESl 
Bi2Se3 has also been studied in a magnetic field, revealing 
a magnetoelectric coupling^^ and an oscillatory angular 
dependence of the magnetoresistance.l^ Sample growth 
dynamics of Bi2Te3 by means of MBE were monitored by 
reflection high-energy electron diffraction (RHEED).^^ 
Doping Bi2Te3 with Mn has been shown to result in the 
onset of ferromagnetism,®'^ while doping Bi2Se3 with Cu 



leads to superconductivity.'^^'^ 

Considerable efforts have been devoted to scanning 
tunneling microscopy (STM) and spectroscopy (STS), 
which enable the study of quasiparticle scattering. Scat- 
tering off surface defects, in which the initial state inter- 
feres with the final scattered state, results in a standing- 
wave interference pattern with a spatial modulation de- 
termined by the momentum transfer during scattering. 
These manifest themselves as oscillations of the local 
density of states in real space, which have been seen 
in several materials with topologically protected surface 
states. Roushan et alP^ have imaged the surface states 
of Bio.92Sbo.08 using STM and spin-resolved ARPES 
demonstrating that, as expected theoretically, in the ab- 
sence of spin-independent (time-reversal invariant) dis- 
order, scattering between states of opposite momenta 
is strongly suppressed. Gomes et alW^ investigated the 
(111) surface of Sb, which displays a topological metal 
phase, and found a similar suppression of backscatter- 
ing. Zhang et al^ and Alpichshev et al^have observed 
the suppression of backscattering in Bi2Te3. Theories 
have recently been formulated to describe the quasipar- 
ticle interference seen in STM experiments, accounting 
for the absence of backscattering.EfflESl Related develop- 
ments have seen the imaging of a surface bound state us- 
ing scanning tunneling microscopy.'^ These experimental 
studies provide evidence of time-reversal symmetry pro- 
tection of the chiral surface states. 

Despite the success of photoemission and scanning tun- 
neling spectroscopy in identifying chiral surface states, 
signatures of the surface Dirac cone have not yet been 
observed in transport. Simply put, the band structure of 
topological insulators can be visualized as a band insula- 
tor with a Dirac cone within the bulk gap, and to access 
this cone one needs to ensure the chemical potential lies 
below the bottom of the conduction band. Given that 
the static dielectric constants of materials under investi- 
gation are extremely large, approximately 100 in Bi2Se3 
and 200 in Bi2Te3, gating in order to bring the chemical 
potential down is challenging. Therefore, in all materials 
studied to date, residual conduction from the bulk exists 
due to unintentional dopiM. One potential way forward 
was opened by Hor et al.^ who demonstrated that Ca 
doping (w 1% substituting for Bi) brings the Fermi en- 
ergy into the valence band, followed by the demonstra- 
tion by Hsieh et al^ that Ca doping can be used to tune 
the Fermi energy so that it lies in the gap. Checkelsky et 
aLl^ performed transport experiments on Ca-doped sam- 
ples noting an increase in the resistivity, but concluded 
that surface state conduction alone could not be respon- 
sible for the smallness of the resistivities observed at low 
temperature. Eto et al^ reported Shubnikov-de Haas 
(SdH) oscillations due only to the 3D Fermi surface in 
Bi2Se3 with an electron density of 5 x lO^^cm"^, while 
Analytis et al^ used ARPES and SdH oscillations to 
probe Bi2Se3. The two methods were found to agree for 
bulk carriers, but while ARPES identified Dirac cone-like 
surface states, even at number densities of the order of 



lO^^cm"^ SdH oscillations identified only a bulk Fermi 
surface, and the surface states were not seen in trans- 
port. A series of recent transport experiments claim to 
have observed separate signatures of bulk and surface 
conduction. Checkelsky et al^ performed weak local- 
ization and universal conductance fluctuation measure- 
ments, claiming to have observed surface conduction, al- 
though the chemical potential still lay in the conduction 
band. Steinberg et aZ.!^ reported signatures of ambipolar 
transport analogous to graphene. However, in this work 
also, two contributions to transport were found, the bulk 
contribution was subtracted so as to obtain the surface 
contribution, and the surface effect was observed only for 
topgate voltage sweeps, not back gate sweeps. Butch et 
aLl^ carried out SdH, ordinary Hall effect and reflectivity 
measurements and concluded that the observed signals 
could be ascribed solely to bulk states. No experimen- 
tal group has reported SdH oscillations due to the surface 
states, which are expected at a very different frequency^^ 
The above presentation makes it plain that, despite the 
presence of the Dirac cone, currently no experimental 
group has produced a true topological insulator. Conse- 
quently, all existing topological insulator/metal systems 
are in practice either heavily bulk-doped materials or thin 
films of a few monolayers^ (thus, by definition, not 3D.) 

In an intriguing twist, two recent works h ave re ported 
STM imaging of the Landau levels of BiaSea.^SiJSaJ Cheng 
et alW^ claim to have observed Landau levels n analogous 
to graphene, with energies ex ynB but the expected spin 
and valley degeneracy of one, yet puzzling features such 
as the absence of Landau levels below the Fermi energy 
and enhanced oscillations near the conduction band edge 
have not been resolved, and some arbitrariness is involved 
in extracting the Landau level filling factor v. Hanaguri 
et al.^"^ report similar findings in the same material, with 
analogous features in the Landau level spectrum. Follow- 
ing the pattern set by the imaging of the surface Dirac 
cone, the Landau levels observed in STM studies have not 
been seen in transport. These last works illustrate both 
the enormous potential of the field and the challenges to 
be overcome experimentally. 

The focus of experimental research has shifted to the 
separation of bulk conduction from surface conduction, 
hence it is necessary to know what type of contributions 
to expect from the surface conduction. The only trans- 
port theories we are aware of have been Ref. [531 which 
focused on magnetotransport in the impurity- free limit, 
and Ref. HH which focused on spin-charge coupling. A 
comprehensive theory of transport in topological insula- 
tors has not been formulated to date. This article seeks to 
remedy the situation. We begin with the quantum Liou- 
ville equation, deriving a kinetic equation for the density 
matrix of topological insulators including the full scatter- 
ing term to linear order in the impurity density. Peculiar 
properties of topological insulators, such as the absence 
of backscattering (which leads to Klein tunneling and 
other phenomena), are built into our theory. We deter- 
mine the polarization function of topological insulators. 



By solving the kinetic equation we determine all contri- 
butions to the conductivity to orders one and zero in 
the transport scattering time r. We find that the scalar 
part of the Hamiltonian makes an important contribu- 
tion to the conductivity, which has a different carrier 
density dependence than the spin-dependent part, and 
may be a factor in distinguishing the surface conduction 
from the bulk conduction. We consider charged impu- 
rity as well as short-range surface roughness scattering 
and determine the density-dependence of the conductiv- 
ity when either of these is dominant. Most importantly, 
we find that the electrical current generates a spin po- 
larization in the in-plane direction perpendicular to the 
field, providing a unique signature of surface transport. 
Moreover, the ratio of the single-particle scattering time 
to the transport scattering time provides important in- 
formation on the dominant angles in scattering, which 
in turn provide information on the dominant scatter- 
ing mechanisms. At zero temperature, which constitutes 
the focus of our work, we expect a competition between 
charged impurity scattering and short-range scattering. 
For each of these two scattering mechanisms the ratio of 
the single-particle scattering time to the transport scat- 
tering time has a qualitatively different dependence on 
the Wigner-Seitz radius rg, which parametrizes the rela- 
tive strength of the kinetic energy and electron-electron 
interactions. We study these aspects of the transport 
problem in detail, covering the possible situations that 
can be encountered experimentally with currently avail- 
able samples and technology. Finally, we consider the 
effect on the topological surface states of adding metal 
contacts, as would be the case in experiment. 

Our work shows that the suppression of backscattering 
seen in experimental work indicates protection against 
localization resulting from backscattering, rather than 
protection against resistive scattering in general. The 
conductivity due to the surface states is determined by 
scattering of carriers by disorder and defects, the density 
of which is in turn determined by sample quality. Strong 
scattering will lead to low mobility, and the mobility is 
not a topologically protected quantity. To see surface 
transport one therefore requires very clean samples. 

Several features set topological insulators apart from 
graphene: the twofold valley degeneracy of graphene is 
not present in topological insulators, while the Hamilto- 
nian is a function of the real spin, rather than pseudospin. 
This implies that spin dynamics will be qualitatively very 
different from graphene. Furthermore, the Dirac Hamil- 
tonian is in fact a Rashba Hamiltonian expressed in terms 
of a rotated spin. Despite the apparent similarities, the 
study of topological insulators is thus not a simple matter 
of translating results known from graphene. Due to the 
dominant spin-orbit interaction it is also very different 
from ordinary spin-orbit semiconductors. 

The outline of this paper is as follows. In Sec. [H] we 
introduce the effective Hamiltonian for the surface states 
of topological insulators and briefly discuss its structure. 



Next, in Sec. HI we derive a kinetic equation for the 



spin density matrix of topological insulators, including 
the full scattering term in the presence of an arbitrary 
elastic scattering potential to linear order in the impu- 
rity density. This equation is solved in Sec. |IV[ yielding 
the contributions to orders one and zero in the trans- 
port scattering time t, and used to determine the charge 
current. Sec. |V] discusses our results in the context of 
experiments, focusing on the electron density and impu- 
rity density dependence of the conductivity, ways to de- 
termine the dominant angles in impurity scattering, and 
the effect of the contacts on the surface states. We end 
with a summary and conclusions. 



Transport experiments on topological insulators seek 
to distinguish the conduction due to the surface states 
from that due to the bulk states. As mentioned in the in- 
troduction, it has proven challenging to lower the chemi- 
cal potential beneath the bottom of the conduction band, 
so that the materials studied at present are not strictly 
speaking insulators. It is necessary for ep to be below 
the bulk conduction band, so that one can be certain 
that there is only surface conduction, and at the same 
time the requirement that Sft/H ;:§> 1 is necessary for 
the kinetic equation formalism to be applicable. These 
assumptions will be made in our work. 



II. EFFECTIVE HAMILTONIAN FOR 
TOPOLOGICAL INSULATORS 

The effective Hamiltonian describing the surface states 
of topological insulators can be written in the form 



H, 



Ok 



Dk^ + A(T-k. 



(1) 



It is understood that k — {k^, ky). Near fc = 0, when the 
spectrum is accurately described by a Dirac cone, the 
scalar term is overwhelmed by the spin-dependent term 
and can be safely neglected, the remaining Hamiltonian 
being similar to that of graphene, with the exception that 
(T represents the true (rotated) electron spin. On the 
other hand, at finite doping and at usual electron densi- 
ties of 10^^ — lO^^cm^^ the scalar term is no longer a small 
perturbation, and should be taken into account. The 
scalar and spin-dependent terms in the Hamiltonian are 
of the same magnitude when Dkp = A, corresponding to 
kp « 10^ m~^ in Bi2Se3, that is, a density of lO^^cm^^, 
which is higher than the realistic densities in transport 
experiments. Consequently, in this article the terms ex D 
will be treated as a perturbation. The eigenenergies are 
denoted by e\ = Dk^ + \Ak, with A = ±. We do not 
take into account small anisotropy terms in the scalar 
part of the Hamiltonian, such as those discussed in Ref. 
1501 Terms cubic in k in the spin-orbit interaction are 
in principle also present, in turn making the energy dis- 
persion anisotropic, but they are much smaller than the 
fc-linear and quadratic terms in Eq. M and are not ex- 
pected to contribute significantly to quantities discussed 
in this work. 

There are contributions to the charge current from the 
scalar and spin parts of the density matrix. The current 
operator is given by 



{2Dk + A(T) 



(2) 



Unlike graphene, where the effective Dirac Hamiltonian 
is valid even at large electron densities and the current 
operator is independent of k for all ranges of wave vector 
relevant to transport, in topological insulators the scalar 
term in the Hamiltonian quadratic in k also contributes 
to the current, and will yield a different number-density 
dependence. 



III. KINETIC EQUATION 

We wish to derive a kinetic equation for the system 
driven by an electric field in the presence of random, 
uncorrelated impurities. To this end we begin with the 
quantum Liouville equation for the density operator p, 



dp 
di 



[i/,p]-0. 



(3) 



where the full Hamiltonian H = Hii+ Hp + f7, the band 
Hamiltonian Hq is defined in Eq. (hj) , Hp ~ eE ■ f repre- 
sents the interaction with external fields, r is the position 
operator, and U is the impurity potential. We consider 
a set of time-independent states {|fes)}, where k indexes 
the wave vector and s the spin. The matrix elements of 
p are p^k' = p^fc/ = {ks\p\k's'), with the understand- 
ing that pkk' is a matrix in spin space. The terms Hq 
and He are diagonal in wave vector but off-diagonal in 
spin, while for elastic scattering in the first Born approx- 
imation U^l., = Ukk'^ss'- The impurities are assumed 
uncorrelated and the average of {ks\U\k' s') {k' s'\U\ks) 
over impurity configurations is {ni\Ukk>\'^Sss')/V, where 
Ui is the impurity density, V the crystal volume and Ukk' 
the matrix element of the potential of a single impurity. 
The density matrix pkk' is written as pkk' — fk Skk' + 
gkk', where fk is diagonal in wave vector (i.e. fkk' oc 
^kk') while gkk' is off-diagonal in wave vector (i.e. k ^ k' 
always in gkk'-) The quantity of interest in determining 
the charge current is fk, since the current operator is 
diagonal in wave vector. We therefore derive an effec- 
tive equation for this quantity by first breaking down the 
quantum Liouville equation into 



dfk 
dt 

dgkk' 
dt 



[HQk,fk] 

[H,g]kk' 



-l[H^Jk]-l[U,g]kkl4^) 



[Uj + g] 



kk' 



(4b) 



and solving Eq. ( 4b ) to first order in U 
gkk'^-T / dt'e-'"'"' Uj(t-t') 



iHt'/h 



kk' 



(5) 



The details of this derivation have been included in Ap- 
pendix [A] We are interested in variations which are slow 
on the scale of the momentum relaxation time, conse- 
quently we do not take into account memory effects and 
f{t - t') « f{t), so that fk obeys 



dfk , « r„ r 
,, + fc [J^Ok,Jk 

at n 


] + Jifk) = 


-I [HE 


Jk], 


(6a) 


Jifk) ^ ^y^ dt' 


U^^-^m'/n 


uJit) 


^tHt'/h 


i^-^) 



The integral in Eq. (6b), which represents the scattering 
term, is performed by inserting a regularizing factor e"''* 
and letting 77 -^ in the end. We consider a potential 
|C^fcfc'| oc 1, which does not have spin dependence, and 
carry out an average over impurity configurations, keep- 
ing the terms to linear order in the impurity density rii. 
The scattering term has the form 



J{fk) = ^ lim 



J2nr 



\Ukk' 



r°° , r 

■J dt'e-^'''{e-^"^'''/^{f,-f,] 



jHkt'/h 



(7) 



For this term to be cast in a useful form one needs to 
complete the time integral and analyze the spin structure 
of the scattering operator J{fk), which is done below. 



A. Scattering term 

The fe-diagonal part of the density matrix /^ is a 2 x 2 
Hermitian matrix, which is decomposed into a scalar part 



and a spin-dependent part. We write fk = n^Tl + S^, 
where n^ represents the scalar part and 1 is the iden- 
tity matrix in two dimensions. The component Sk is 
itself a 2 X 2 Hermitian matrix, which represents the 
spin-dependent part of the density matrix and is written 
purely in terms of the Pauli a matrices. Thus, the matrix 
Sk can be expressed a.s Sk — -^ Sk ■ cr = -^ SkiCi- With 
this notation, the scattering term is in turn decomposed 
into scalar and spin-dependent parts 



JU) 



d^k' 



h^ J (27r)2 
d^k' 



2h^ J (27r)2 



\Ukk'\^{nk-nk^) / dt'e-''*'e-^'='*'/''e'^'=*7/^ + /,.c. 



\Ukk' 



{Sk - Sk') ■ / dt' e-''*'e-*^'='*'/'"'cr e'""''/'' + h.c. 
Jo 



(8) 



We assume electron doping, which implies that the chemical potential is in the conduction band £+. For the scalar 
part Uk and spin-dependent part Sk of the density matrix the above results in 



-iHyt' /h^iHkt' /h _|_ y^^ 



■nh 



{l + k-k')[5{e\-e+)+5{e'_-e^)] 



+ !I^^.(^k + y)[S{e'^-e+)-5{e'_-eJ] 
e-*^fc'*7 V e'^"''/" + h.c. -^ ^ [(t(1 - k ■ k') + [k ■ (t) k' + k [k' ■ a)] [S{e'+ - e+) + 6{e'_ - e_ 

+ ^ik + k')[S{e'^-e+)-Sie'_-e^)], 



(9) 



where fe is a unit vector in the direction of fc. We emphasize that there is no backscattering, so the physics of topological 
insulators is built into the kinetic equation from the beginning. Collecting all terms, we obtain the scattering term 



as J if) = J(n) + J(S'), with 



J(n) = 


nui f d k' 


2h J (2^)2 


+ 


TTUi f d?k' 


2h J (2^)2 


J{S) = 


TTUi f d?k' 


4/1 J (2^)2 




■Kn, f £k' 



Ah J (27r) 



lUkk'^rik ~ nk'){l + k-k') [S{e'^ - e+) + d{e'^ - e_)] 

|C/fefefK - nk')(T ■ {k + k') [5{e\ - e+) - <5(e'_ - e_)] 

\tJkk'?{Sk - Sk') ■ [(t{1 -k-k') + {k-a)k' + k {W ■ rr)] [5{e\ - e+) + 5[e'_ - e_ 

|C7fcfcf (5fc - Sk') ■ {k + k') [5{e'+ - 6+) - 6{e'_ - e_)]. 



(10) 



The details of this calculation are contained in AppendixlAJ The energy (5-functions 5(e'j_ — e±) = [\ / {A±2Dk)] S{k'—k). 
The (5-functions of e_ are needed in the scattering term m order to ensure agreement with Boltzniann transport. This 
fact is already seen for a Dirac cone dispersion in graphene where the expression for the conductivity found in Ref. [SSI 
using the density-matrix formalism agrees with the Boltzmann transport formula of Ref. [86] (the definition of r differs 
by a factor of two in these two references.) The necessity of keeping the e_ terms is a result of Zitterbewegung: the two 
branches A = ±1 are mixed, thus scattering of an electron requires conservation of e-|_ as well as e_. In the equivalent 
kinetic equation for graphene there is no coupling of the scalar and spin distributions because of the particle-hole 
symmetry inherent in the Dirac Hamiltonian. In topological insulators the Dk^ term breaks particle-hole symmetry 
and gives rise to a small coupling of the scalar and spin distributions. Later below we quantify this coupling by an 
effective time td, which — >■ cx) as 13 — >■ 0. 

We find a series of scattering terms coupling the scalar and spin distributions. Let 7 represent the angle between k 
and fe', and the unit vector 6 in polar coordinates represent the direction perpendicular to fe, that is the tangential 
direction. The unit vectors fe, A;', 6 and 6' are related by the transformations 



k' = k cos 7 + sin 7 
0' = —ksinj + 0COS7 



k = k' cos J — 6' sin 7. 



(11) 



We decompose the matrix Sk = Sk\\ + Sk± and write those two parts in turn as Sk\\ = (1/2) Sfc|| cr^w and Ski_ = 
(1/2) Sfcj^cTfcj^. The small 5^11 and Sfc_L a-re scalars and are given by Sk\\ = S^ ■ k and Sk± = 3^0, and similarly 
ffcll — cr ■ k and ak± = cr ■ 6. We introduce projection operators P„, P\\ and P±^ onto the scalar part, ak\\ and crfcj_ 
respectively. We will need the following projections of the scattering term acting on the spin-dependent part of the 
density matrix 



P\\J{Sk\\) = 
Pj{Sk\,)^ 
P\\J{Sk±) = 



krii 


1 


1 


IGh-K 


XA + 2Dk) ' [A- 


-2Dk) 


kfii 


1 


1 


lefiTT 


{A + 2Dk) ' [A- 


-2Dk) 


kn.i 


1 


1 


le/iTT 


{A + 2Dk) ' {A- 


-2Dk) 



dO' \Ukk'\^ (sfcii - Sfc'ji)(l + cos 7) crfcii 



dO' \Ukk'\'^ {sk\\ -Sfc'||)sin7crfe_L 



d9'\Ukk'f {sk± + Sk'±) sin j(Tk\\ 



(12) 



and PA_J{Ski_) will be deferred to a later time. Furthermore, the scattering terms coupling the scalar and spin 
distributions are 



J, 



n^>-S 



(nk) 



Js^n{Sk) — 



7mj 

~2h 
4h 



(Pk' 
d^k' 

(2^ 



\Ukk'\^{nk - rifc')[fTfc|[(l + COS7) + crfe_Lsin7] [5{e'j^ 



-)-5{e'_-e^)] 



(13) 



\Ukk'?[{sk\\ - Sfc'l|)(l + cos7) + {ski_ -hSfc'_L)sin7] [^(e^ - e+) - (5(e'_ -£-)]. 



Notice the factors of (1 -I- cos 7) which prohibit backscat- 
tering (and lead to Klein tunneling.) 



B. Polarization function and effective scattering 
potential 

The above form of the scattering term is valid for any 
scattering potential, as long as scattering is elastic. To 



make our analysis more specific we consider a screened 
Coulomb potential, and evaluate the screening function 
in the random phase approximation. In this approxima- 
tion the polarization function is obtained by summing 
the lowest bubble diagram, and takes the form^^ 



fcAA' 



/ofcA — fok'X' 

\- £fcA — £fc'A' - 



irj 



1 + AA' cos 7 



(14) 

where /ofcA = foiskx) is the equilibrium Fermi distribu- 
tion function. The static dielectric function, of interest 
to us in this work, can be written as e{q) ~ 1 + v{q) n('?), 
where v{q) = e^ / (2eoerq) ■ To determine n(g), we assume 
T = and use the Dirac cone approximation as in Ref. 
[87l This approximation is justified by the fact that we 
are working in a regime in which T/Tp <C 1, with Tp 
the Fermi temperature. The results obtained are in fact 
correct to linear order in _D. At T = for charged impu- 
rity scattering the long- wavelength limit of the dielectric 
function ia^ 



e{q) ^ 1 






(15) 



yielding the Thomas-Fermi wave vector as kxp — 
e^kp/iineoerA). 

As a result, in topological insulators the matrix ele- 
ment Ukk' of a screened Coulomb potential between plane 
waves is given by 



Ukk' 
\Ukk'\' = ^^ { r, . ,: , , 1 (16) 



Ze' 




1 


2eoe, 


\k- 


k'\ + kTF 


Z^e^ 


( 


1 


^44 


\\k 
w 


-k'\+kTP 



2fei? / 



where Z is the ionic charge (which we will assume to be 
Z — \ below in this work) and k'pp is the Thomas-Fermi 
wave vector. The parallel projection of the scattering 
term becomes 



P,lJ(5,l) = 



kuiW 
16fnr 

krii 
8h^ 



1 



1 



{A + 2Dk) {A-2Dk) 
1 1 



de' 



1 -I- cos 7 

■ 7 I kTF,-' 



{s\\ -S|i)'^fe|l 



C(7) = 



{A + 2Dk) {A-2Dk) 

2 7 



2 ' 2kp I 

d0'c(7)(sii -4)^^11 



(17) 



cos 



We expand C(7) ^-nd ^k\\(P^ in Fourier harmonics 

C(7)=EC™e^"^ 

in 



AmO 



(18) 



P\\J{S\\ 



1 



1 



kui 

^ [{A + 2Dk) ' {A-2Dk) 

r 



X!(^0-Cm)s||me" 



O'fclh 



and to obtain the expansion of Sk\\{d') replace 6 — ^ 9'. 
The absence of backscattering, characteristic of topolog- 
ical insulators, is contained in the function C(7)- Equa- 
tion (17) shows that this function contains a factor of 
(1 -|- cos 7) which suppresses scattering for 7 = tt. 



IV. SOLUTION OF THE KINETIC EQUATION 

The Hamiltonian Hp describing the interaction with 
the electric field is given by Hp = eE ■ r, with r the 



position operator. It then follows straightforwardly that 
the kinetic equation takes the form 



dfk 
dt 



[HkJk] + J{fk)='D- 



(19) 



dpoh 
dk 



The driving term in the kinetic equation isT) — ^ 
where pofc represents the equilibrium density matrix, 

given by pofc = (/ofc+ + /ofc-) l + (l/2) (/ofc+-/ofc-) cr-fc- 
The first term represents the charge density, while the 
second represents the spin density. We decompose the 
driving term into a scalar part 2?„ and a spin-dependent 



part Vs, as 
2?„ 



eE-k ,a/( 



0fe4 



dfi 



Ok- 



dk 



dk 



1 fcE-k dfok^ 



a/c 



Ofc- 



dk 



(20) 



1 







eE-e 



(/( 



Ofc4 



/ofc- ) • 



2~ ~ hk 

The latter is further decomposed into a part paraUcl to 
the Hamiltonian Vu and a part perpendicular to it Vj^ 



Vn 



1 eE-k .dfi 



Ok-i 



dk 



dfok- 
dk 



1 



lo-fell 



J-k\\ O'fell 






1 



/Ofe- ) CTfci = X C'fei CTfci • 



(21) 

Evidently, since wc consider electron doping, /ofc- = 0. 
Nevertheless we retain /ofc- for completeness, since the 
term (/ofc+ — /ofc-) will be shown to give rise to a con- 
tribution to the conductivity singular at the origin. 

We solve the kinetic equation. The first step is to di- 
vide it into three equations: one for the scalar part, one 
for the part parallel to the Hamiltonian and one for the 
part orthogonal to the Hamiltonian. 

duk 



dSk 



+ PnJifk) 

P\\J{fk) 



= T^n, 



= v« 



(22a) 
(22b) 



dS, 



fc_L 



dt 



dt 

[Hk,Ski]+Pj{fk) = V^. (22c) 



We search for the solution as an expansion in the small 
parameter h/{£FT), where t cx n^ , to be defined below, 
represents the transport relaxation time. We label the or- 
ders in the expansion by superscripts with e.g. SjT^ rep- 
resenting Sk evaluated to order m. From these equations 
it is evident that the expansions of Uk and Sk\\ begin at 
order Spr/h, in other words [h/{eFT)Y~^' , whereas Sk± 
begins at order [h/{eFT)Y^\ in other words the leading 
term in it is independent of epT /h. Therefore, in deter- 
mining the leading terms in Uk and Sfcy, the coupling to 
Sfc_L can be neglected. We are then left with two coupled 
equations for Uk and Sk\\, which are written as 

Jn^n{nk) + Js^n{Sk) = 'D-n 



1 dSk\\ 

2 "dT 



c^fell 



dt 



P\\J: 



s^s{S\\) + P\\Jn^s{nk) = 2 4|| o-fc|| 



(23) 

We must carry out the projections in order to determine 
the scattering terms above. In the equations involving P|| 
we have eliminated (Jk\\ , since the equations for Sk have 
been cast in terms of equations for Sk- To first order in 
D, the required scattering terms can be written as 

J{n)^^-^jde'C{l){nk-nk.) 
T /ON kuiWDk f ,„, ^, ■, , ■, 

Js^n[Sk) -^ , t,A2 / "-^ C(7) (Sfell - Sfc'll) 



P\\Jn^sink) 



kuiWDk 
27rM2 



dd' (ij) {uk -Uk') (7k\\ 



The scattering terms have exactly the same structure. 
We proceed to solve the equations for nk and Sk\\- We 
note that the RHS on both equations in this case has 
only Fourier component 1, so we can make the Ansatz 
Uk = ni e~*^ -I- n_i e'^ and Sk\\ = syi e~*^ -I- S||_i e*^. 
Observing in addition that I?„ = dk\\ allows us to write 
the coupled equations in the steady state as 



2nfc 

T 

2nfc 

TD 



Sk\\ 
TD 



2 ' 



(25) 



where t, the transport relaxation time, is given by the 
expression 



1 kmW 



2hA 



(Co-Ci), 



(26) 



and we have also introduced 



1 kn,WDk 



td 



hA^ 



(Co - Ci), 



(27) 



as a measure of the coupling of the scalar and spin distri- 
butions by scattering. The equations are trivial to solve 
and we obtain in the steady state 



nk = ^il + T/2TD) 
Sku = ^ (1 + 2T/Tn) 



(28) 



The equation is valid to first-order in t/te,. Note that, 
as D — > 0, r£) — >■ oo and we recover the results for the 
bare Dirac cone. 

Next we examine the term Sl,j_. The kinetic equation 
for 4°| is 



jc(0) 

-^ + I Wk, 4"|] =d^- P^Anit'] - PJ[Sit\ 

(29) 
The projections of the scattering terms do not contribute 
to the conductivity. The bare source term d±_ is 



«_L — — ^fTJ— UOfe+ - JOfc- 



hk 



(30) 



The result for Sl,l_ is 



S; 



(0) 
k± 



eE-e 

2hk 



(/( 



0/c+ 



lOk- 



dt e ''* ( ak± cos 

I 



2Akt 



(T ■ k X 9 sin 



2Akt 



(31) 



The part ex cr ■ k x 9 averages to zero over directions 
in momentum space and is not considered any further. 
Integrating the cosine over time gives 



(^0) _ eE-9 
^k± - 2hk ^ ^ 



V 






O-fcJ, 



(32) 



A. Current operator and conductivity 

The current density operator is decomposed into a par- 
allel and a perpendicular part, the electric field is as- 
sumed II X and we are looking for the diagonal conduc- 
tivity 



IMPLICATIONS FOR TRANSPORT 
EXPERIMENTS 



In this section, as in the calculation of screening above, 
it is assumed the system has a low enough density that 
the Dirac-cone approximation is applicable, which for 
Bi2Se3 implies kp < 10®m~^. Below we use the termi- 
nology high density as meaning that the density is high 
enough that the Dk^ term in the Hamiltonian has a no- 
ticeable effect (meaning n w lO^^cm^'^ and higher), with 
low density reserved for situations in which this term 
is negligible. Nevertheless we retain the scattering time 
T£) , which allows us to examine the correction to the low- 
density expressions when these terms become noticeable. 



ji{k) = jt^n + - jj,|| Cr|| + - j^.± cr_L- 



(33) 



The three contributions to the current operator, repre- 
senting the scalar part, the part parallel to the Hamil- 
tonian and the part orthogonal to the Hamiltonian, are 
given explicitly by 



Jx,n 



Jx,\\ 



Jx,± 



2eDk 
2eA 



cost 



cost 



(34) 



2eA 



suit 



Using this form for the current operator it is trivial to 
show that the conductivity linear in r, keeping only terms 
linear in £), is given by 



.Boltz 



e 



AkpT 
Ah 



2t 
td 



Dkl 



Ah 



(35) 



We emphas ize ag ain that this result tends to the result 
for graphen d^^ l ^^ l when 13 — >■ (note that the definition of 
T in Refs. 18515^ differs by a factor of 2.) The contribution 
to the conductivity independent of r is 



,0_L 



7re 



lim 



1 



1 + p/{tJ.+ h^/2) I ^ ^P(jj.-hLo/2) 



(36) 
The latter term however is zero for finite /i, which is the 
case we are considering. Thus, the leading order term 
in the conductivity is oc r, and in the limit epr/h ^ 1 
there is no term of order t*'. 



A. Density dependence 



The conductivity is a function of two main parameters 
accessible experimentally: the carrier number density 
and the impurity density/scattering time. The depen- 
dence of the conductivity on the carrier number density 
arises through its direct dependence on kp and through 
its dependence on n through r. In terms of the number 
density the Fermi wave vector is given by kp = Awn. The 
transport time is given by Eq. ( 26 1 and W is defined in 
Eq. (16). The number-density dependence of r depends 



on the dominant form of scattering and whether the num- 
ber density is high or low as defined above. For charged 
impurity scattering W ex A:~^, so that r ex fc and the 
two terms in the conductivity are ex n and n'^". At the 
same time, in a two-dimensional system surface rough- 
ness gives rise to short-range scattering as discussed in 
Ref. [HHl For short range scattering of this type W is 
a constant and r ex fc, with the two terms in the con- 
ductivity being a constant and n^/^. These results are 
summarized in Table IT] At high density the term ex Dk 
becomes noticeable through the scattering time td. Since 
Tp) has an extra factor of kp, retaining terms to first or- 
der in D gives us an additional term ex Dn^'^, resulting 
in a further contribution to th e co nductivity of the form 
ex Drt'/^ , as is manifest in Eq. (35 1 . Taking this fact into 
account, in order to keep the expressions simple, in Ta- 
ble H] we have listed only the number density dependence 
explicitly, replacing the constants of proportionality by 
generic constants. 



10 



TABLE I: Carrier density dependence of the conductiv- 
ity. High-density and low-density regimes for scattering by 
screened charged impurities and for short-range impurities. 
The numbers ai — 04 represent constants. 

Screened charges Short — range 



Low density 
High density 






■il2 



03 
tta + 04?! 



1/2 



TABLE H: Ratio of the transport and Bloch scattering times 
T JTQ for Ts ^ 4, Ts ^ 4 and different scattering mechanisms. 
61 — feg are constants. 



Screened charges Short — range 



<4 
»4 






65 -l-&6''slnra 
h^ -I- hars 



Ukk' = U{q) as 



fcpn,; 



AnhA 



kprii 



tq AttKA 



dj 



dj 



Ujq) 

4q) 



(1 — cos^ 7) 



(1 -f- COS 7). 



(38) 



We recall that for elastic scattering of carriers on the 
Fermi surface |fc — fc'| = q — >■ 2fci?sin7/2. We consider 
first charged impurity scattering, denoted by the sub- 
script c. The two scattering times can be expressed in 
terms of a prcfactor with the units of frequency multi- 
plying a dimcnsionless angular integral, thus 1/tqc = 
{rl/Toc)lQc and l/r^ = {r^/Toc)Itc, where the Wigner- 
Seitz radius rg — e^ /(2TTeQerA), the quantity 



B. Scattering times 

Two scattering times commonly encountered in mo- 
mentum relaxation: the transport relaxation time r and 
the quantum lifetime tq . The former represents momen- 
tum relaxation in transport and is weighted by a factor of 
(1 — COS7) to reflect the fact that forward scattering does 
not alter the charge current. This is the time measured 
in transport experiments. The latter represents the time 
taken to scatter from one Bloch state into another, thus 
all scattering amplitudes are weighted equally. This time 
determines observables such as the Dingle temperature 
and the broadening of Shubnikov-de Haas oscillations. 
As a result of the different weighting of the scattering 
amplitudes in r and tq, the ratio t/tq characterizes the 
scattering by providing information on which angles are 
dominant.!^ We determine the ratio of these two scat- 
tering times in topological insulators. First, we provide 
expressions for r and tq in the presence of charged im- 
purity scattering. To begin with, we write the quantum 
lifetime tq and the transport lifetime r in a slightly dif- 
ferent way from the way they were written in Sec. 
(withD = 0, and d = 0): 



IV 



1 



kpTiiW 
T ^ AttHA 

1 kpriiW 



dlCil) (1 ~ cos 7) 



■^Q 



AttHA 



27r 



(37) 



d-fCil), 



where tq is discussed at length in Ref. [HH We can also 
express them in terms of the original scattering potential 



1 

TOc 



Afij 



(39) 



and the angular integrals It and Iq 



Itc{x) 
Iqc{x) 



dj 



dj 



1 — cos 7 
(sin ^ + x^ 

1 + cos 7 



(40) 



(sin ^ + xy 



Even before evaluating the integrals explicitly it will be 
observed that the transport time t is not sensitive to ei- 
ther small or large angle scattering, while the quantum 
lifetime tq is insensitive to large angle scattering. The 
qualitatively different dependence on the scattering an- 
gles makes plain the fact that the ratio t/tq will have 
a qualitatively different behavior as a function of r^ for 
charged impurity scattering and for short-range scatter- 
ing. It must be emphasized also that, despite the super- 
ficial similarity to graphene, the behavior of this ratio as 
a function of r^ for real materials is also quite different, 
due to the absence of the spin and valley degeneracy, and 
the fact that Ts for a topological insulator differs from its 
value in graphene (for the same parameters) by an overall 
factor of 2. 

Substituting for the screening function we can deter- 
mine the expressions for these times as a function of r^ 
in the limits of small and large r^. The results are qual- 
itatively the same as in Ref. |92l and are summarized in 
Table [n] For charged impurity scattering, we find for the 
integral Iqc 
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Iqc{x) 



^ (2 — a;^)arccoth-\/l — x'^ 

.-2" + ^^ (l-X^)3/2 



27r,x = 1 
4 8x 



1 — a;^)arccothW J— ^] 



,x <1 



-27T + 



X Vx2^^ 



ar 



rccot V a;2 — 1 — arctan i / ) , a; > 1 



2: + 1 



X- 1 



(41) 



while Ifr can be written as 



Itc{x) = 2{tt + Vlx ~ %-KX^) + 
32 - IOtt, a; = 1 

2(7r + Vlx - eTra:^) - 



16a;(3x2 - 2) 



16a;(3x2 - 2) 



ar 



ccothv 1 — x"^ — arctanh 



1 +.T 
1-X 



a; < 1 



arccot v x'^ — 1 — arctan \ I — ' — ] ,x > 1 



x + 1 
x-1 



(42) 



Setting X ^ rg/A we find that the ratio of the transport 
and quantum hfetimes for charged impurity scattering in 
the limit r^ ^ 4 can be expressed as 



^Q 



8 
7rr., 



16 In) 



(43) 



For short-range scatterers characterized by a potential 
vo6{r), the transport and quantum scattering times are 
given by the following expressions 



n tqs ' 4 ' 

TQ5 Tqs 4:' 

where the quantity tqs is defined by 



1 niVgy/n 



Tqs 2Ay/TTh 



(44) 



and the integrals appearing in Eqs. (44 1 are 



hs{x) = i d-i 







(1 — cos^ 7) 



Iqs{x)= / ^71:^1^^^) (I + COS7). 







(45) Evaluating the integrals explicitly yields 



(46) 



Iqs{x) = tt + 127ra; — 6t 



16 — 57r, a; = 1 



TT + 127ra; — 67ra; 



8x(3a:^ — 2)(arccoth-\/l — x^ — arccothw jr^) 



-,a; < 1 



8a;(3a:^ — 2)(arccoth\/a:^ — 1 — arccothw Iq:^) 



VT^^ 



-,x > 1, 



(47) 



and 



Itsix) 



n 16a; 
2 ^ ^ 
1 
6 



(256-8l7r),a; = l 



vr 16x 



Gttx^ + AOx^ - 207rx^ 
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16a:;'^(5a;^ — 4)(arccoth\/l — x^ — arccothw jr^) 



,a; < 1 



(48) 



16a;'^(52;^ — 4)(arccot\/a;2 — 1 — arccotw |— j) 



V^2T7i; 



,a; > 1. 



Setting again x = rg/A the ratio of the transport and 
quantum hfetimes in the hmit Ts ^ 4 for short-range 
impurities takes the form 



250 






4r^hir, 2(17 - 181n2)r^ 



(49) 



In the range covered here this ratio does not exceed 2 for 
short-range scatterers. We would also like to point out 
that Tg ex. 1/A for the surface of Bi2Te3 is approximately 
twice that in graphene, i.e. r^ « 0.14, which is <C 4.1221 



C. Numerical results 

We conclude our discussion with numerical results on 
the conductivity and the ratio of the two scattering times. 
From the transport scattering time calculated in the pre- 
vious section we can calculate the conductivity using Eq. 
(35) for d = 0. In Fig. [T] we have plotted the conductiv- 
ity versus the ratio n/rii of the carrier number density to 
the impurity number density for several values of d, the 
distance of the impurities from the surface of the topo- 
logical insulator. For d = 0, when the impurities are 
located right on the surface, the conductivity is a linear 
function of n/rii. If d is large, meaning kpd > 1 then 
large angle scattering is strongly suppressed. This leads 
to a large increase in t/tq, since the Fourier transform 
of the scattering potential is suppressed by e~'^'^. 

In Fig. [2] we have plotted the ratio of the two scatter- 
ing times as a function of the carrier number density for 
charged impurity scattering for different values of the dis- 
tance d of the impurities from the surface. For d — this 
ratio is not a function of n. It should also be noted that, 
for short range scatterers, the ratio ts/tqs is a constant 
as a function of n, which we find to be 1.46. 

We note that the analytical results determined in the 
previous sections for d = agree with the numerical re- 
sults shown in Figs. 1 and 2. Nevertheless, for d 7^ 0, 
it is not possible to compare directly the numerical re- 
sults with the analytical results, the latter being valid for 
d — 0, i.e. the impurities located on the surface. 

Ref. |80] measured a mobility of the order of 1 m^/(Vs) 
at a number density of the order of « lO^^cm"^. For 
example, for one sample the surface carrier density is es- 
timated at ss 7 X lO^^cm"^, while the effective scattering 




n/ii: 



1 



FIG. 1: Calculated conductivity limited by screened charged 
impurities as a function of the surface carrier density for dif- 
ferent impurity locations d = 0, 5, lOA. In this figure we 
use the following parameters: an impurity density rii = 10^^ 
cm~'^, j4 = 4.1 eVA, which corresponds to the Fermi veloc- 
ity vf = 6.2 X 10^ cm/s, and the static dielectric constant 
er = 100. 



rate is 0.8ps, which is obtained from transport. For such 
a sample epT /h ~ 35, placing it within the range of ap- 
plicability of our theory. However the experiment is not 
sensing surface state transport so a direct comparison is 
not feasible at present. 

Our theory has focused on scattering due to charged 
impurities and interface roughness, whereas in a realistic 
sample other scattering mechanisms may exist, primarily 
phonons. Phonons can be included easily in our theory as 
an additional contribution to the scattering term. How- 
ever in this work we have considered zero temperature, 
therefore there is no phonon absorption and there is only 
phonon emission, which usually has a negligible contri- 
bution to the conductivity. 

We note that our theoretical results would be useful in 
discerning transport by the surface topological 2D states 
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FIG. 2: The ratio of the transport relaxation time to the 
quantum life time for screened charged impurities as a func- 
tion of the carrier density for different impurity locations. We 
use the same parameters as Fig. [l] The ratio for the short- 
range impurity is independent of carrier density and we find 
ns/TQs = 1.46. 



only after a real topological insulator material, i.e. a sys- 
tem which is a bulk insulator, becomes available. All 
currently existing TI materials are effectively bulk met- 
als because of their large unintentional doping. Such 
systems are, by definition, unsuitable for studying sur- 
face transport properties since bulk transport completely 
overwhelms any surface transport signature. Discussing 
TI surface transport in such bulk-doped TI materials is 
not particularly meaningful since it must necessarily in- 
volve complex data fitting and assumptions in order to 
distinguish bulk versus surface transport contributions. 
Although such analyses are being carried out in bulk 
doped TI materials, in order to understand the transport 
behavior (and our theoretical results should certainly 
help such analyses), we believe that any real progress will 
come only when surface TI transport can be carried out 
unambiguously, without any complications arising from 
the (more dominant) bulk transport channel. We believe 
that our theoretical results presented in this paper would 
become important in the context of such surface TI stud- 
ies in the laboratory. 

What we have shown in this paper is that the so-called 
topological protection of TI surface transport is a pro- 
tection only against localization by back scattering, not 
a protection against impurity or defect scattering. The 
presence of impurities and defects will certainly lead to 
scattering of the TI surface carriers and the surface 2D 
conductivity will be strongly affected by such scatter- 
ing. If such scattering is strong (which will be a non- 



universal feature of the sample quality) , then the actual 
surface resistivity will be very high and the associated 
mobility very low. There is no guarantee or protection 
against low carrier mobility in the TI surface states what- 
soever, and unless one has very clean surfaces, there is 
little hope of studying surface transport in 2D topologi- 
cal states, notwithstanding their observation in beautiful 
band structure measurements through ARPES or STM 
experiments. 



D. Minimum conductivity at the Dirac point 

Thus far we have concentrated on the conductivity Uxx 
at carrier densities such that epT /h ^ 1. This is the limit 
that is immediately relevant to experiment given the high 
unintentional doping in topological insulators. At such 
densities the conductivity is given by a^°^^^ from Eq. 
(35), while the term a^^ from Eq. (36) does not con- 



tribute. At present it is not possible to tune the number 
density continuously through the Dirac point as is done 
in graphene. Nevertheless, given the similarity of the 
Hamiltonians of graphene and topological insulators, we 
expect a similar behavior at the Dirac point. We predict 
that, when such tuning does become possible, a mini- 
mum conductivity plateau around the Dirac point will 
be observed, analogously to graphene. 

The term cr^x i^ expected to provide one contribu- 
tion to the minimum conductivity, which can also be 
obtained using a field-theoretical method in the absence 
of disorder. ^^ This method can in fact yield several val- 
ues, depending on the choice of regularization proce- 
dure, which are close in magnitude. However, as n — >■ 
the condition spT/h ^ 1 is inevitably violated, and a 
strong rcnormalization of a^^ is expected due to charged 



impurities.' 



Quite generally, charged impurities give 



rise to an inhomogeneous Coulomb potential, which is 
screened by both electrons and holes. The net effects of 
this potential are an inhomogeneity in the carrier den- 
sity itself and a shift in the Dirac point as a function of 
position. At high densities the spatial fluctuations in the 
carrier density are of secondary importance, and do not 
modify the linear dependence of the conductivity on n, 
yet as the chemical potential approaches the Dirac point, 
where the average carrier density (n) = 0, these fluctu- 
ations play the dominant role in conduction.®^ At low 
densities, the carriers cluster into puddles of electrons 
and holes, and a residual density of carriers is always 
present, making it impossible to reach the Dirac point 
experimentally.'^ Consequently, although Eq. ( 36 ) ap- 



pears to give a universal value for the minimum conduc- 
tivity, the renormalization due to the presence of electron 
and hole puddles displays a strong sample dependence: 
in graphene it res ults in an enhancement of a^^ by a 
factor of 2-20. '^^'^ In addition, exchange and correlation 
effects make a signiflcant contribution to the minimum 
conductivity.'^ 

An accurate determination of this enhancement for 
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topological insulators requires detailed knowledge of the 
impurity density distribution, which can only be deter- 
mined experimentally or modeled numerically by means 
of e.g. an effective medium theory:- ^ However, a self con- 
sistent transport theory provides a physically transparent 
way to identify the approximate minimum conductivity. 
The magnitude of voltage fluctuations can be calculated 
in the random phase approximation, and the result used 
to determine the residual density and the critical num- 
ber density at which the transition occurs to the regime 
of electron and hole puddles, where the carrier density 
will be highly inhomogeneous. Both the residual density 
and the critical density are proportional to rii by a factor 
of order unity,!^ and as a first approximation one may 
take both of them as ~ n,. In view of these observations, 
the minimum conductivity plateau will be seen approx- 
imately at <J^°''*^, in which with the carrier density n is 
replaced by rii. Using the expressions found above for 
CT^^'*^ and r, this yields for topological insulators 



(50) 



At carrier densities approaching zero the ZJ-dependent 
term in the Hamiltonian will not contribute. Therefore, 
when experiment will be able to tune the carrier density 
to zero we expect a minimum conductivity of the order of 
cr™" to be observed. We note that, with r^ much smaller 
than in graphene due to the large dielectric constant, the 
minimum conductivity may be substantially larger.EB 



E. Current-induced spin polarization 

The charge current is proportional to the spin operator, 
as can be seen from Eq. ^ . Therefore a nonzero steady- 
state surface charge current automatically translates into 
a nonzero steady-state surface spin density. The spin 
density can easily be found by simply multiplying the 
charge current by ti^/{—2eA), yielding 



in 



AkpT 
4h 



2t 
td 



Dkl 



Ah 



(51) 



Since cr in the Hamiltonian represents the rotated spin, 
this corresponds to a net density of Sy. The effective 
Hamiltonian of Eq. (fTl) is obtained from the original 
Rashba Hamiltonian by replacing Uy — > Ux and Ux — >■ 
—(Ty For a sample in which the impurities are located 
on the surface, with n/rii = 0.5 and Ex = 25000^/771, the 
spin density is « 5 x 10"'^* spins /m? (where spin = h), 
which corresponds to approximately 10~'* spins per unit 
cell area. This number, although small, can be detected 
experimentally using a surface spin probe such as Kerr 
rotation. It is also a conservative estimate: for very clean 
samples, having either a smaller impurity density or im- 
purities located further away from the surface, this num- 
ber can reach much higher values. Current-induced spin 
polarization is a signature of two-dimensional transport 



in topological insulators since there is no spin polariza- 
tion from the bulk. 



VI. EFFECT OF METALLIC CONTACTS 

A question of major practical importance concerns the 
stability of the surface states in the presence of metallic 
contacts. Transport measurements require the addition 
of metallic contacts on the surface of the topological in- 
sulator. Since the properties of the surface states depend 
crucially on the boundary conditions, the natural ques- 
tion is how are these properties affected by the metallic 
contacts. To address this question, we consider the case 
of large area contacts and perform numerical calculations 
for the minimal tight-binding model studied in Ref. [Ml 
A topological insulator in contact with a metal along a 
planar interface is described by the Hamiltonian 



^tot — H- 



TI 



H 



M 



Ht 



(52) 



The topological insulator term i?Ti is given by the tight 
binding model on a diamond lattice with spin-orbit 
interactions^ 



Hr 



xr 



-t 



J 






((u)) 



..(di^xd^JcLc 



■jo-': 



(53) 



where the first term is the nearest neighbor hopping on 
the diamond lattice and the second term connects second 
order neighbors with a spin and direction dependent am- 
plitude. The direction dependence is given by the bond 
vectors d^,' traversed between sites i and j. The metallic 



term in Eq. (52 I is given by 



H, 



M 



{ij) ,o- 



(54) 



The model is defined on a simple hexagonal lattice with 
a lattice constant for the triangular basis that ensures 
simple interface matching conditions between the metal 
and the (1, 1, 1) surface of the TI. The tunneling term in 
Eq. ([52l) is 






(55) 



were t is the tunneling matrix element that characterizes 
the transparency of the interface between the topological 
insulator and the metal. 

For a free surface, the topological insulator spectrum 
is characterized by a bulk gap occupied in the vicinity 
of the M-point by an anisotropic Dirac cone associated 
with gapless surface states."^ These surface states are 
exponentially localized near the boundary, with a char- 
acteristic length scale of a few lattice spacings. When the 
TI surface is in contact with a metal, the surface states 
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FIG. 3: (Color online) (a) Amplitude of a metallic state hy- 
bridized with a topological insulator surface state as function 
of distance from the interface (in units of interlayer spacing). 

(b) Spectrum of a topological insulator in contact with a metal 
(black) . The states of a topological insulator with a free sur- 
face are also shown (purple - bulk states, red - surface states). 

(c) Density plot of the total amplitude within a thin region of 
topological insulator in the vicinity of the interface (red/dark 
grey area in panel a). Note that, instead of a well-defined 
surface mode, one has a diffuse distribution of states with 
boundary contributions. The dispersion of the maxima of 
this distribution is represented by green points in panel (b). 



penetrate inside the metal and hybridize with the metal- 
lic states. A typical hybridized state is shown in Fig. [3^. 
Note that the amplitude of this state near the boundary 
of the topological insulator is reduced by a factor 1/Lz 
with respect to the amplitude of a pure surface state with 
the same energy, where Lz is the width of the metal 
in the direction perpendicular to the interface. However, 
the local density of states near the boundary is not re- 
duced, as the number of hybridized states also scales with 



Lz ■ The spectrum of a topological insulator in contact 
with a metal is shown in Fig. |3]d. Also, we calculated 
for each state the total amplitude within a thin layer of 
topological insulator in the vicinity of the boundary. The 
result is shown in Fig. |3J2. Instead of the sharply defined 
Dirac cone that characterizes the free surface, one has 
a diffuse distribution of states with boundary contribu- 
tions. 

The properties of the interface amplitude distribution 
shown in Fig. |3j:, i.e., its width and the dispersion of 
its maximum, depend on the strength of the coupling 
between the metal and the topological insulator. If the 
distribution is sharp enough and the dispersion does not 
deviate significantly from the Dirac cone, our transport 
analysis should be applied using bare parameters, i.e., 
parameters characterizing the spectrum of a free surface. 
Nonetheless, significant deviations from the free surface 
dispersion are, in principle, possible. For example, the 
location of the Dirac point in our model calculation is 
fixed by symmetry. However, in real systems the en- 
ergy of the Dirac point can be easily modified by chang- 
ing the boundary conditions. Consequently, the effective 
parameters -Dgs and A^s entering transport coefficients 
may differ significantly from the corresponding param- 
eters extracted from ARPES measurements. To clarify 
these issues, as well as to determine the dependence of 
these effects on the size of the contacts, further studies 
are necessary. 



VII. SUMMARY 

We have determined the charge conductivity of the 
surface states of three-dimensional topological insulators 
basing our theory on the quantum Liouville equation and 
accounting fully for the absence of backscattering. We 
identify two contributions to the conductivity with dif- 
ferent dependencies on the carrier number density. At 
low densities the conductivity is ex n if charged impurity 
scattering is dominant and is independent of n if short- 
range scatterers such as surface roughness are dominant. 
As the density increases an additional term is present in 
the conductivity, which is oc n?'"^ for charged impurity 
scattering and ex n^/^ for short-range scattering. The 
ratio of the transport and quantum lifetimes in the re- 
gion Ts <C 4, which we believe to be the most realistic 
experimentally, is of the form hirj^ — 62 Inrc, for charged 
impurity scattering and of the form 65 + b^rglnrg for 
surface roughness scattering. The analysis we have pre- 
sented should enable one to obtain a complete character- 
ization of all the features of topological insulator samples 
accessible in transport experiments, provided £pT /h ^ 1 
and the chemical potential is below the bottom of the 
conduction band. 

The current form of the theory presented in this arti- 
cle gives results equivalent to those of Boltzmann trans- 
port theory. Nevertheless, this theory accounts explicitly 
for interband coherence effects such as Zitterbewegung, 
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which should become important near the Dirac point in 
ballistic samples. We expect the advantages of this for- 
mulation to become manifest in the presence of a mag- 
netic field. A magnetic field induces additional coherence 
between the two bands, influences the spin dynamics, and 
gives rise to a qualitatively different scattering term. The 
flexibility of the theory outlined in this work enables us 
to study a variety of multiband effects beyond Boltzmann 
transport, which we reserve for a future publication. 
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Appendix A: Derivations of certain expressions appearing in the scattering term 

Firstly, Eq. (l5| is derived as follows: 

dgkk' , * r/y M * rr> f , -1 

-^ + -^[H,9U, = --[UJ + gU,, 



(Al) 



To first order in U we do not retain the g-dependent term on the RHS. We also drop the k and k' indices for simplicity, 
writing the equation as an operator equation 



Now let g = e^^^^^'^gie^^'^^^, a transformation equivalent to switching to the interaction picture 



(A2) 



dgi 



'— e 

dt h 



gi 



iHt/h 



[Uj]e 



-i.Ht/h 



— CXD 

oo 



dt"e~''"'"/\U, f{t - t")y"'"/^ 



(A3) 



Projecting onto k and k' we obtain Eq. (Is]). 

We consider in more detail the derivation of the contribution oc (fc + k') in the scattering term. The scattering 
term has the general form (with summation implied over the index i) 



Jif) 



K^ 



J w? '^'"'' '' ^"^ " "''^ C '^^' ^'''^'^"'"' ''^"'' + ^-^ 



+ ^J§^ \Ukk'\'{Sk^ - ^.'0 J^^ dt' e-'"'e-^'='*'a, e^"-'' + h.c. 



(A4) 



In the time evolution operator one substitutes the general expression for the band Hamiltonian Hof, = Dk^ + (1/2) cr ■ 
Jlfe, where fl^ — 2Ak. In what follows we shall shorten the time evolution operators by omitting the factors of h. 
These will be restored in the final expressions. The time evolution operator can be written as 



r^fci 



Hut 



-iHkt —iDk''tl '^'fe'- ■ A • ^^k 

e '° = e '^'^ ' I cos — I cr • rife sm — 



(A5) 



We recall the identity 



(cr • rjfc) (cr • tlk' ) — tlk- tlk' +i(T -tlky. ^h 



(A6) 
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This allows us to express the products of two time evolution operators as 

ttlkt ttl^it ^ l(UK -t^K ;tr g (,Qg Z<T-JlfcCOS Sm hiCT-Jlfc/Sm COS 

^2 2 22 2 2 

+ (S7fc/ • l7fc + zcr • Slfc X Slfc'j sm — — sm ——] 

e ^^^Ve *= = e '- •• [crcos ^^ cos — « (cr • fife) <t cos — — sm — V irr (ct ■ Hfc-) sm — — cos — - 

, f f\ ^ t A ^ ■ ^k't . ^kt-, 
+ {cT ■ ilk) a- [a- ■ ilk') sm — — sm — -J. 

The time integrals of these products of exponentials and trigonometric functions give 

^.(Dk'^^Dk^)t ^^3 ^ eos ^ ^ ^ [5(e^ .\) + 5{e^ - e'_) + 5{e^ - e'_) + 5{e^ 4)] 

^.iUk^-Dk^n eos %^ sin ^ ^ J [S(e'^ - ._) + S{e'_ - e'_) - S{e'^ .+) - 6{e'_ .+ )] 



(A7) 



(A8) 



The Hermitian conjugates are trivial. For the scalar part of the density matrix, letting O^ = k the above results in 



+ ^^.(fc + fc')[<5(4-e+)-'5(e--e-)]. 



(A9) 



For the spin-dependent part of the density matrix, 

e-'"^'a- e'""'' ^ ^ [cr + (o- • fc) cr (cr • k')] [5{e'^ - e+) + 5{e'_ - e_)] 
+ ^ [(^ . fc) ^ + <T (<T . fe')] [<5(e'+ - ^+) - '^(e'^ - £_)] 
e-*-f^'='*'cr e^^-"*' + e^*^'"*'^- 6*-^'='*' ^ — {[2cr + (cr • fc) cr (cr • fc') + (cr • k') a {a ■ k)] [5(e'^ - e+) + &{e'_ - e_)] 

+ — [{cT-k)a + (T{a- k') + {a-k')a + (T{(T- k)] [6{e'+ - e+) - ,5(e'_ - e_)]}. 

(AlO) 
The product of three sigma matrices simplifies as follows 

(cr • fc) cr (cr • k') + (cr • k') cr (cr • fc) = 2(fc ■ a) k' - 2{k ■k')(T + 2fc (fc' • cr), (All) 

yielding finally 

e-'"'''*' (T e'""'' + e-^^'-'V 6*'^'='*' ^ ^ [cr(l - k ■ k') + {k ■ a) k' + k (fc' • cr)] [J(e'^ - e+) + 5(e'_ - e_)] 

(A12) 
+ ^(fc + fe')[<5(6;-6+)-<5(e'_-6_)]. 



18 



The scattering term takes the form J(/) ~ J[n) + J{S) 
(Pk' 



Jin) = ^ y (^ It/fefcf (nfe - nfc,)(l + fc • fc') [S{e'^ ~ e+) + J(e'_ - e_)] 



+ ^ / ^ It/fefcf (nfc - nfc> . (fe + fc') [5(6; - 6+) - <5(e'_ - 6_ 



2?i y (27r)2 



J(5) 



^ / (S l^'^'^'l'(^'^ - ^'^') • ['^(1 -k-k') + {k-cT)k' + k {k ■ a)] [5{e'+ - e+) + J(e'_ - e_)] 
/ ^ l^'^'^'I'l^'^ - Sw) ■ (k + k') W^ 6+) - 5{e'_ 6_)]. 



(A13) 



TTTl. 

^ Ui J (27r) 

Note the factor of 2 difference between J{n) and J{S), which arises solely from our definition iJofc = Dk'^ + {l/2) erClk- 
Note also that there are no second-order terms in the Pauli matrices ai. From the algebra of these matrices we know 
that 



(JiCTj — dij + eijkCTk 

[(Ti,(Tj] — 2eijk<Tk 

{ai,aj} — 2Sij. 



(AM) 



in other words the second-order terms can give zero, the identity matrix [leading to the terms oc (Sk — SkA ■ {k + k')], 
or one Pauli matrix. 
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